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We consider the Ising phase of the antiferromagnetic XXZ Heisenberg chain on a finite-size lattice 
with N sites. We compute the large N behavior of the spin stiffness, obtaining the correlation length 
£. We use our results to discuss the scaling behavior of metal-insulator transitions in ID systems, 
taking into account the mapping between the XXZ Heisenberg chain and the spinless fermion model 
and known results for the Hubbard model. We study the scaling properties of both the Hubbard 
model and the XXZ Heisenberg chain, by solving numerically the Bethe Ansatz equations. We find 
that for some range of values of £/N the scaling behavior may be observed for the Hubbard model 
but not for the XXZ Heisenserg chain. We show how £ can be obtained from the scaling properties 
of the spin stiffness for small system sizes. This method can be applied to models having not an 
exact solution, illuminating their transport properties. 



PACS numbers: 71.27+a, 71.30.+h, 72.15.Nj 

I. INTRODUCTION 

Since Bethe's solution 1 of the one-dimensional (ID) 
isotropic Heisenberg model, that the Bethe Ansatz (BA) 
method for solving and understanding ID many-body in- 
tegrable systems has been a valuable tool in condensed 
matter. 

Many physical quantities can be computed from the 
BA equations alone. Among these are thermodynamic 
quantities, which can be computed using the thermody- 
namic Bethe ansatz method. 2 Nevertheless, the BA so- 
lution does not provide a way of computing all physical 
quantities of interest. For example, non-diagonal correla- 
tion functions, such as the spectral function, important 
to understand photoemission experiments, or the opti- 
cal conductivity, important to understand the response 
of the system to an electromagnetic field have to deter- 
mined using other methods. 3-5 

The study of transport properties in integrable sys- 
tems is one of the most active fields of research in in- 
tegrable models in these days. The two most studied 
models have been the Hubbard and the spinless fermion 
(or the spin 1/2 Heisenberg chain) models. The under- 
standing of the transport properties of these systems is 
related to the calculation of the optical conductivity (or 
spin conductivity 6 ). This quantity is defined as 

a(ui) = 2ttD c S(u) + a r (ui) , (1) 

where D c is the charge stiffness and a r (u>) is the regular 



part of the conductivity. This last quantity is associ- 
ated with absorption of radiation at finite u> values. The 
charge stiffness is associated with the existence of a metal 
insulator transition in the system, driven by the inter- 
action (a Mott-Hubbard transition). If D c is finite the 
system presents infinite conductivity (free acceleration), 
but if D c — the system is an insulator. This simple 
picture holds at zero temperature. 

W. Kohn showed that D c can be computed from the 
ground state energy, Eq, of the system, imposing twisted 
boundary conditions, as 7 
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where 4> is a twisting angle (or, in the ID flux 
piercing the ring). 

At finite temperature, the picture above is slightly 
changed but Eq. (1) still holds. 8 The first question one 
can raise is whether at finite temperature the system still 
has a finite value of D c if D c is finite at zero tempera- 
ture. The answer to this question was first conjectured 
by Zotos and Prelovsek using numerical methods 9 , and 
put on firm grounds by Zotos 10 in the context of the XXZ 
Heisenberg model. These authors have shown that inte- 
grable systems present ballistic transport of charge and 
spin at finite temperature if the charge (or spin) stiffness 
is finite at zero temperature. 

These theoretical questions about the effect of inte- 
grability on transport properties triggered experimental 
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research on transport of energy and spin in quasi ID 
materials, such as SrCuC>2 and Sr2Cu03. The results 
are controversial. The experimental study of heat trans- 
port by spin excitation of SrCuC>2 and Sr 2 Cu03 materi- 
als revealed an anomalous enhanced thermal conductiv- 
ity for temperatures above 40 K, 11 and the authors con- 
cluded for ballistic transport of energy by the spinons. 
The MNR study, 12 for zero momentum transfer, in the 
Sr2Cu03 material was able to probe the uniform spin sus- 
ceptibility separately from the staggered ±ir/a (a is the 
lattice constant) mode. The authors concluded for ballis- 
tic spin transport in the T = limit and for diffusion- like 
transport at finite temperatures (even for T <C J, where 
J is the antiferromagnetic coupling). Both from an ex- 
perimental and theoretical 13 point of view, the study of 
transport properties in ladder spin systems is open to re- 
search, and the questions raised in this work should also 
be pursued in those systems. 

Although the theoretical picture for the spin trans- 
port in the XXZ Hciscnbcrg model is clear if we con- 
sider the system in the parameters region correspond- 
ing to the XY universality class [A < 1, see Eq.(3)], a 
complete understanding of the spin stiffness (or charge 
stiffness, in terms of the spinless fermion model) in the 
parameters region corresponding to the Ising university 
class (A > 1) is lacking. In the XY universality class 
region the zero temperature spin stiffness is given by 
D c = J7rsin(^i)/(4^(7r — /j,)), with 1 < [i < n, reach- 
ing the value J/4 at the isotropic point (p, = ir). 14 At 
finite temperature, and for 1 < [i < n, the spin stiffness 
was computed by Zotos, 10 although there is some contro- 
versy about its temperature dependence. 15 In the Ising 
universality class region the zero temperature spin stiff- 
ness is zero, in the thermodynamic limit, and its value 
for finite temperatures has not been computed yet but it 
was conjectured it should also be zero. 9 ' 16 We stress here 
that the above conjecture is by no means obvious. For 
example, in the Hubbard model at half-filling, the zero 
temperature charge stiffness is zero, but it was found to 
be finite at finite temperatures, 18 although this result 
has been criticized by computing the curvature of levels 
explicitly in the U ^> t limit. 19 

Also the scaling properties of D c in the region A > I 
have not been studied as were in the A < 1 region, 17 and 
for the Hubbard model. 20 

In the study of the scaling properties of the Mott- 
Hubbard insulator the existence of a scaling function 
Y(£/N) for the charge stiffness was shown, 20 where £ is 
the correlation length and N is the number of sites, and 
it was conjectured that Y(£/N) should be a universal 
function characterizing other metal-insulator transitions 
in ID systems and, in particular, it should characterize 
that transition in the spinless fermion model. 

The spinless fermion model can be mapped into 
the spin 1/2 Heiscnberg chain by a Jordan- Wigncr 
transformation. 21 Therefore the transition between a spin 
conductor (spin currents propagating ballisticaly) and a 
spin insulator (existence of spin diffusion), occurring at 



the isotropic point, and characterized by the spin stiff- 
ness value can be mapped into a metal-insulator transi- 
tion, driven by the Coulomb interaction among spinless 
fermions at adjacent sites. In this work we consider the 
spin language, but our results are equally valid for spin- 
less fermions. 

This work is organized as follows. In the section II 
we derive the asymptotic (large N) behavior for the spin 
stiffness and obtain the correlation length £ charactering 
the insulating state. In section III we present a numerical 
study of the spin stiffness of the XXZ Heisenberg chain 
and the Hubbard model. In section IV we summarize our 
results. 



II. ASYMPTOTIC BEHAVIOR 



The Hamiltonian of XXZ model reads 



n = -jfl^ s r + i e ~ iHN + H - c -) + ASf sf +1 ], (3) 



where N is the number of sites, , Sf , Sf are spin 1/2 
operators at site I, and is a flux piercing the ID ring. 
The connection with spinless fermion model is made by 
means of the Jordan- Wigner transformation. In terms of 
spinless fermion the Hamiltonian (3) can be written as 
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where the spinless fermion operators c\ , Cj obey the usual 
anti-commutation relations, means summation over 
nearest neighbors, n, = c\ci is the usual local number op- 
erator, ^ is the hopping integral (normally called t but 
here called ^ to stress the important exact correspon- 
dence between this model and the XXZ Heisenberg spin 
1/2 chain), and V = J A is the nearest-neighbor Coulomb 
repulsion. 

The Hamiltonian (3) can be solved by the BA, 2,22,23 
and the energy eigenvalues are given, for A = cosh(7) 
( 7 > 0), by 
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e = e q -jY^ 
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where E = JAY/4, N a (a =|, J.) stands for the num- 
ber of up and down spins, and the numbers dis- 
cussed below. In the case ATj = N± the system has an 
energy gap relatively to the ground state. In terms of 
spinless fermions the system is an insulator (in the ther- 
modynamic limit). The energy gap is given by 



00 (— I)™ 
J sinh 7 > — — 



cosh(n7) 



(6) 



2 



The numbers Xj obey Bethe ansatz equations (with 
twisted boundary conditions) 
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N6i {xj , 7) = 2-irIj +4> + J2 d 2 ( x j -xul). (7) 

1=1 

where 



n (x, 7) = 2 tan 1 



tan(7x/2) 



tanh(n7/2) 

The ground state quantum number configuration is: 
Ni-1 Ni-3 Ni-1 



(8) 



(9) 



In the absence of magnetic flux, the ground state den- 
sity is given by 



1_ A y, cos(ri7a:) \ 
\Q \ ^ cosh(n 7 ) J 



2Q 
K 



■dJ^ 



2ttQ V 



(10) 



where K — K{k) is the complete elliptic integral 24 of 
the first kind, dn(y, k) is a Jacobian elliptic function of 
modulus fc 24 and we have 
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The ground state density (10) can be obtained from the 
momentum density p(x) =am{Kx/Q,k) as dp(x)/dx = 
2np(x), and where am(y, fc) is the amplitude function. 

In order to obtain the charge-stiffness large- N asymp- 
totic behavior we write the energy (6) as 
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J2S(x- Xj ). (12) 



Using the properties of the delta function and the Poisson 
summation formula, Eq. (12) can be written as 



E = E Q - JN sinh 2 7 



p N {x,cj>) dx 



cosh 7 — cos(7a;) 

00 

exp{im[Np N (x) - 4>]} , (13) 
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where pN(xj) = (2nlj + <j>)/N. For large N the integral 
in Eq. (13) is dominated by the m = 0, ±1 terms. As a 
consequence the stiffness (2) is a sum of two terms and 
reads 



D c = D^ m=±1) + Di m=0) 



(14) 



where 
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The first term dominates Eq. (14) for large values of 
A, since the second one is, in this case, exponentially 
small. Both Eq. (15) and Eq. (16) are evaluated using 
the steepest descent method. The final asymptotic result 
for D r can be cast in the form 



D C {N) ~ D(-y)y/Ne 



(17) 



where £(7) = ip(Q + 1) — inN/2 is the correlation length 
given by 
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and D("f) is given by 
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The correlation length (18) agrees with the value ob- 
tained previously using a different method, 25 and is rep- 
resented in Fig. 1. 
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FIG. 1. Representation of the correlation length as func- 
tion of A. It is clear that £ diverges as we approach the 
isotropic point A = 1. 
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The result (18) is quite similar to that obtained for the 
Hubbard model by Stafford, Millis, and Shastry. 20 



III. NUMERICAL RESULTS 

At T = the XXZ Heisenberg model for A > 1 
presents non-ideal spin transport in the thermodynamic 
limit, because the spin stiffness is zero. In terms of spin- 
less fermions one has an insulating behavior. For a finite- 
size system D c is finite even for large values of N and for 
moderate values of A, as can be seen in Fig. 2. In this 
figure we see that even for values of N as large as 1000 
sites there are values of A for which D c is of the order 
of its value at the isotropic point (represented by a filled 
circle). This shows that the behavior of D c , when A > 1, 
in a finite-size system, is dominated by the value of the 
D c for A = 1. The value of D c at the isotropic point is 
by no means trivial, and the presence of this point has 
consequences for the scaling behavior of D c in the Ising 
region as we show below. We believe this difference is due 
to the presence of logarithmic corrections to the scaling 
function by marginal interactions at the isotropic point. 



0.3 




1 2 3 

A 



FIG. 2. Representation of the stiffness, D c , as function of 
A for different system sizes. The value of D c at the isotropic 
point (A = 1) in the thermodynamic limit is represented by a 
filled circle. The represented values of A satisfy the condition 
A > cosh(0.4). 

The above discussion stresses that for mesoscopic sys- 
tems the separation between an insulator and a metal is 
not clear. Only when the correlation length is smaller 
than the system size a clear separation between an insu- 
lator and a metal exist. The correlation length (18) is 
independent of the system size, and the dependence of 
D c on £ is represented in Fig. 3. This figure makes the 
above statement clear, for example, for values of £ bigger 
than 10 2 , D c is finite even for N = 1002, whereas for 
values of £ smaller than 10 2 D c is zero for N = 1002 but 
not for the other (smaller) values of N. 
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FIG. 3. Representation of the stiffness, D c , as function of 
£ for different system sizes. The values of £ correspond to 
values of A satisfying the condition A > cosh(0.4). 

In Ref. 20 the charge stiffness of the Hubbard model, 
at half-filling, was shown to present a scaling behavior of 
the form Y(£/N), for all values of the Hubbard interac- 
tion U, where £ is the correlation length and N is the 
system size. In Figure 4 we show data from a numeri- 
cal solution of the BA equations for the Hubbard model, 
which agrees with the results of Ref. 20. Looking more 
carefully to Figure 4 we see there is scaling behavior for 
all values of N as we get closer to the gapless phase U = 
(£ — ► oo). On the other hand, as U increases, £ — > 0, and 
the scaling behavior can only be verified for larger system 
sizes. 
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FIG. 4. Representation of the stiffness, ttD c , as function of 
£/N for the Hubbard model, for system sizes N =24, 32, 42, 
62, 82, 100, 302 and 502. The correlation length £ is given by 
Eq. (3.2) of Ref. 20. 

Based on the scaling behavior of the Hubbard model, 
a conjecture that the same scaling behavior would ap- 
ply to other ID systems (in particular for the spinless 



fermion system) presenting metal-insulator transitions 
was drawn. 

Following these ideas we have first derived an asymp- 
totic expression for the stiffness of the XXZ Heisenberg 
model (or spinless fcrmions). This allowed us to deter- 
mine the correlation length £ for this model. We have 
then solved numerically the BA equations for systems of 
different sizes and tried to see if some type of scaling was 
possible. From Figure 5 we see that the representation of 
D c as function of £/7V gives curves that do not collapse 
into each other unless N > 202. This larger value of N 
need for scaling to be observed contrast with the case of 
the Hubbard model in Figure 4, although the values of 
£/N are the same in both cases. Also, as we get closer 
and closer to A = 1, where the system presents a gapless 
phase (in Figure 5 this corresponds to £/N — > oo), the 
collapsing of the data takes place for larger and larger val- 
ues of N. This behavior contrasts with that seen above 
for the Hubbard model. 
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FIG. 5. Representation of the stiffness, D c , as function of 
£/N, for the Heisenberg model. The inset shows the same 
but for N =502 and 1002, being clear that the collapsing of 
the data close to A = 1 requires larger values of N. 

Since the scaling as function of £/i\T is not perfect for all 
values of N, let us assume that D c could follow a scaling 
behavior of the form D c = N~ n f({;/N) and let us find the 
value of the exponent r\ that leads to the collapsing of the 
data. We proceed by grouping the several N in groups 
of two as follows {22, 50}, {50, 102}, {502, 1002}, and 
for each set r] is determined by a minimization procedure 
which leads to the best possible scaling. Our findings are 
listed in Table I, and show that the value of the exponent 
77 decreases as the N increases, leading for very large val- 
ues of N, to the predicted scaling behavior D c = f(£/N). 
On the other hand, for small values of N the scaling be- 
havior D c — N~ v f(^/N) works much better as we show 
in Figure 6. The fact that it was not possible to find 
single scaling behavior for for all N is an indication that 
either a different scaling law or corrections containing N 
and associated with the presence of the critical point are 



need. 

Let us now show that scaling relation D c = 
N~ v f(£/N) can be used to determine £ directly from 
the solution of small systems. The scaling relation 

<7W0 = ^/(£) (20) 

has an interesting property which can be used in order 
to determine the correlation length £ in cases where this 
can not be done analytically. Supposing that (20) holds 
for a given quantity, g, and assuming the existence of a 
correlation length depending solely on other parameter 
of the problem, i. e. £ = £(A) (which, in the case of our 
problem, is the anisotropy A), the scaling relation can be 
expressed in a slightly different, but equivalent notation 
as 

gx {N) = g(N,t{\)) = ±f(Qp) . (21) 

We can now consider, as is frequently the case in nu- 
merical calculations, measuring g for various sizes and 
different values of the parameter A. Then, for each pair 
of sets g Xi (N) and g Xj (N) we can write 

9K vm) = m [m) f ' (22) 

or, setting nj = 

g Xj {N)=g Xi {Nr ij )rl . (23) 

This relation should hold for given nj and 77, that can be 
extracted by direct minimization 26 . Proceeding in this 
way for the entire range {Ai,...,Am} in which the nu- 
merical values of g are calculated, the dependence of £ 
in the parameter A can be obtained, except for a global 
multiplicative constant. An example of the typical result 
obtained via this procedure can be seen in Figure 7 (we 
remark that for this case the value of 77 is constant for 
all values of N and A used). Here, numerical results for 
the Heisenberg model with chains of sizes 6 to 24 were 
used to calculate the charge stiffens D C (N,A). We re- 
mark that the values of N used here are easily accessible 
to exact diagonalization. Although the scaling relation 
(20) doesn't seem to apply exactly to this problem, it was 
found, for A not too close to the critical point, that a 
good scaling can be obtained as shown in the left frame 
of Figure 7, allowing the numerical calculation of £ (A), 
as shown in the right frame of Figure 7. As can be seen, 
we were able to calculate the correlation length without 
any further information other than the numerical results 
for D C (N,A). This concept can be usefully applied to 
other cases where the analytical result for £ is unknown, 
at least to give a rough idea of the behavior of the corre- 
lation length or any other quantity involved in the scaling 
relation. 
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FIG. 6. Representation of the stiffness, N V D C , as function 
of £/N, for N = 22 and TV = 50. 
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FIG. 7. On the left panel we present the collapse of the 
numerical data for several distinct A according to the scal- 
ing relation (21). On the right are our numerical results for 
£ (A) obtained by means of the described method (squares), 
together with the exact result (circles) obtained from the 
Bethe ansatz, for comparison. 

IV. SUMMARY 

We have studied the spin stiffness of the XXZ Heisen- 
berg chain for A > 1. We have shown that the in large 
N/£ behavior of the stiffness agrees with the asymptotic 
behavior D(N) = Y(N/£), where Y(x) is a universal 
scaling function, but not for all ranges of iV/£, as in the 
Hubbard model. Contrary to the case of the Hubbard 
model and for a given value of £/iV the system sizes re- 
quired to observe the scaling behavior are about an order 
of magnitude larger for the XXZ Heisenberg chain. More- 
over as we get closer and closer to the isotropic (A = 1) 
point, characterized by £/N — > oo, larger and larger sys- 
tem sizes are required, indicating the scaling behavior 
does not apply in that region. We believe the difference 



between the Hubbard model and the XXZ Heisenberg 
chain is related to the presence of the non-trivial point 
A = 1, for the Heisenberg chain as opposed to the triv- 
ial point U — (independent electrons) for the Hubbard 
model. An interesting question which remains to be an- 
swered concerns the scaling properties and the correlation 
length of the metal-insulator transition in an extended 
Hubbard model, considering next-nearest neighbor inter- 
action, since this model is no longer integrable. For small 
system sizes the above scaling function does not apply 
and a scaling of the form D(N) = N~' n f(N/£) produces 
a much better scaling behavior. Using this scaling func- 
tion the correlation length £ was obtained from the anal- 
ysis of the small systems alone, and shown to agree well 
with the exact value. This method can be used to derive 
the correlation length for systems which do not have an 
exact solution. 
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